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ABSTRACT 

The group theory framework developed by Fukutome for a systematic 
analysis of the various broken symmetry types of Hartree-Fock solutions 
exhibiting spin structures is here extended to the general many-body context using 
spinor-Green function formalism for describing magnetic systems. Consequences 
of this theory are discussed for examining the magnetism of itinerant electrons in 
nanometric systems of current interest as well as bulk systems where a vector spin- 
oq density form is required, by specializing our work to spin-density-functional 
O formalism. We also formulate the linear response theory for such a system and 
^ compare and contrast them with the recent results obtained for localized electron 
On systems. The various phenomenological treatments of itinerant magnetic systems 
are here unified in this group-theoretical description. 

a 

£ 

I. INTRODUCTION: 

a 

O Irreducible representations of the symmetry group of a system which classify 

the eigenstates of its Hamiltonian have been discussed before [1]. However, many 
important physical phenomena such as superconductivity and many types of 
magnetic states of condensed matter display lower symmetry characteristics, and 
are thus broken symmetry states of the Hamiltonian [1]. These phenomena have 
been theoretically understood in terms of approximate solutions which break the 
full symmetry of the Schrodinger equation of the many - body system. For 
example, the celebrated Hartree - Fock (- Bogoliubov) self-consistent field method 
has been used in this way. A clear description of this approach may be found in [1] 
where an enumeration of the various commonly found basic symmetries is also 
given. Fukutome [2] and his coworkers, based only on symmetry arguments, 
systematized the search for possible novel states of the system, which in the past 
depended on intuitive suggestions ( viz. Bardeen, Overhauser, and Landau, for 
example). These authors developed a complete group-theoretical classification and 
characterization of all the possible magnetic symmetry structures arising from only 
the underlying symmetry group consisting of spin rotations (S), time reversal (T), 
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and spatial symmetry group (P) and its subgroups, within the Hartree - Fock 
approximation. Fukutome found that there are only eight subgroups of the 
symmetry group SxT, consisting of spin rotations and time reversal. The 
Hamiltonian of the well-known electron gas which is an important ingredient in the 
density-functional method, is invariant under SxT. It is common knowledge that 
the electron gas exhibits several states of magnetic order when treated in HF and 
related mean field approximation [3]. The eight different structures found 
correspond to those belonging to the eight subgroups. When the spatial symmetry 
group is combined with these, one gets many more varieties of structures. This 
theory has been more recently extended to include mean - field type analysis of 
Hubbard - type models of one- and two - dimensional systems of great current 
interest [4]. In recent years, spin - density - functional method [5] has become the 
method of choice replacing the Hartree - Fock (HF) method both because it 
includes the important correlation effects not included in the HF - scheme in 
describing the theory of magnetic systems and because of the present much 
improved computing techniques. There has been a suggestion to incorporate the 
Fukutome classification to the spin - density - functional formalism [6]. The Green 
function method for describing the many - body systems is more general than the 
above schemes and the broken - symmetry solutions can be incorporated in this 
framework as noted already in [1]. A generalized version of the density - functional 
method which includes spin as well has been reformulated in the Green function 
language [7] . The various broken symmetry states of the HF scheme will be more 
generally expressed here in terms of the symmetry of the Green function under the 
appropriate subgroups for which we give the nomenclature "broken- symmetry- 
adapted Green functions." The purpose of this paper is two fold; one, to 
incorporate the Fukutome classification into the Green function framework and, 
two, to focus attention on magnetic systems in general, by applying this formalism 
to the spin-density-functional method with an emphasis on the currently interesting 
findings of magnetic states in nanometric systems [8]. By this extension, we also 
provide a more complete vector - spin - density - functional (VSDF) theory of 
magnetic systems which goes beyond the currently used schemes for handling 
noncollinear magnetic materials [9, 10]. The recent theoretical work on the 
magnetism of free atomic clusters exhibiting interesting geometric structures 
accompanying very large magnetic moments [11] provides another class of 
problems in magnetism requiring further investigation particularly if the suggestion 
that these clusters can be formed on semiconductor interfaces [12]. In this case the 
possibilities of variety of geometric structures of the cluster formation each with its 
own magnetic features in relation to the substrate would provide some challenge in 
determining the appropriate broken symmetry states of magnetism in such 
environment. The group theoretical enumeration would make this search 
systematic and orderly. 

In Section II, we first give a brief description of the magnetic system in 
terms of the spinor Green function [13] and the associated spinor self - energy 
operator [14]. We then develop the Fukutome classifications of magnetic systems 
associated with this spinor Green function and the related spinor self - energy (or 
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mass -) operator, which is more general than the effective self - consistent potential 
of the HF theory. In Section III, the linear response functions arising out of these 
Green functions [15] are examined to establish the stability aspects of the broken 
symmetry structures. This will be shown to be related to the recent work [16]. In 
Section IV, the implications of the above developments to the VSDF theory are 
spelled out. This is complementary to the work in [17]. In the final Section V, a 
summary and concluding remarks in relation to the proposed work on magnetism 
of nanostructure systems including magnetic atomic clusters are given. 



II. THE BROKEN SYMMETRY ADAPTED GREEN FUNCTIONS - THE 
FUKUTOME CLASSIFICATIONS IN MAGNETIC SYSTEMS 

(a) Hamiltonian: 

The general spin polarized system without spin - orbit interaction (for 
simplicity of presentation at this juncture) is described by the standard system 
Hamiltonian 

H = T e + V u + V ie + V ee . (1) 
Here, T e is the kinetic energy operator of the electrons, V u , V ie , and V ee are 
operators representing the Coulomb interactions between the ions (i) and the 
electrons (e). Introduce y/ s (r), the usual field operator annihilating an electron of 

spin s at position r, and y/l(r), the corresponding creation operator. Then we 
have the definitions of electron density operator 

h{r) = s Q {r) = £ y/l(r)¥s(h (2) 



and the vector spin - density - operator 

hr) = Y,* s M(r)V s (h ( 3 ) 

s,s' 

where the Pauli spin matrix vector is T. It is useful to consider the system as being 
subjected to external, spin-dependent field described by 

Vex = L 1 dr ( rt) yrt ( r) ^ ( r). (4) 

s,s' 

In the following we split w ss < into a scalar part w n = f which acts on the electron 
density given by eq.(2) and a traceless part w s expressed in terms of the Pauli 
matrices which acts on the spin density, eq.(3): 

w s {rt)=i-f{rt). (5) 
We thus rewrite eq.(4) in a physically transparent form 

Vex, = \d?[h{r)w n {rt) + s(r) ■ f(rt)_ 



= H\drs a {r)f a {rt\ 



(6) 



a 

where we have introduced the 4- vector notation 
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(r)-K,/), (tJ-(1,t), and 



t ,r\_ ,r\ (7) 



In the above Greek indices run from to 3. The Hamiltonian (1) is invariant under 
all spin rotations, S, and time-reversal, T. 

(b) Spinor-Green's Function and the Spinor-Self Enegy: 

The one-particle Green's function is written as a matrix in spin indices [13], 

G ss , {rt;¥ f ) = -i(r( ¥s {rt) wl (? f ))}• (8) 

The equation of motion satisfied by this spinor-Green function is often written in 
its most general form as 



d V? 



i +-L-y c (l)-w(l) 
at x 2 m 



G(12)-Jj(3)i:(13)G(32) = 5(12). (9) 



Here the spin-independent classical Coulomb potential Vc arising from the 
electron density n(l) = = -iTrG(l 1 + ) and the nuclei (Vi n) is given by 

V c {\) = \{d2)v{\2)n{2) + V ion {\), (10) 
where the instantaneous Coulomb interaction between electrons is 

v(12) = r - r 8{t x -t 2 ). Also, Tr stands for trace over spin indices only. Also 



r l~ r 2\ 



w(l) = w„(l)+ T-/(l), is the spinor representation of the external field. We have 

used here 1 to stand for the space-time point and the other notations are 

standard usage as in [13], for example. The last term in the left hand side of eq.(9) 
is the contribution due to inter-particle interaction in its most general form and we 
call "exchange-correlation spinor-self-energy" contribution, in anticipation of later 
application to functional formalism. The general expression for this spinor-self- 
energy matrix, Z , is in general a functional of the spinor-Green function, G, and is 
given by 

£(12) = ij d(3)j rf(4)G(13)A„(32;4)W(41) (11) 
with A n (12;3) = -^^ (12) 

is the screened charge-response vertex function, with V tot = V c + w n . W is the 
Coulomb interaction v(12), screened by the dielectric function, e, both of which 
are spin scalars, defined by 

e- 1 ^) d3) 

dw n {2) 

W{\2) = } d(3)£-\l3)v(32). (14) 
Finally we write the Dyson equation (9) in the familiar form 
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G" 1 (12) = 



fA + ZL-y c (i)- w (i) 

2m 



<5(12)-E(12) 

(15) 



- Go 1 - E. 



V c + Z = —, (17) 



In the usual HF approximation, E is the familiar unscreened exchange self-energy 
if we do not include the screening by the dielectric function of the medium: 
E(12) « E# F (12) = /v(12)G //F (12). This exhibits explicit dependence of the self 
energy on the Green function and hence the self-consistency feature of eq.(15). In 
[13], the spinor structure of G was used to construct a self-consistent solution to 
the general HF approximation in the Green function language, from which the 
general Overhauser's spiral spin-density-wave solution was deduced. It was then 
recognized as a "broken symmetry" solution of the HF equation because this 
solution breaks the spin-rotation and time-reversal symmetries of the interacting 
electron gas Hamiltonian! We should also point out that in constructing any 
approximation scheme for the self-energy E, it is useful to have certain 
conservation principles as well as variational character. Baym [14] set up such a 
conserving scheme by introducing a functional <&[G] whose first variational 
derivative with respect to G yields the exact E and constructed an expression for 
the grand potential 

Q[G] = Trtr{\n(-G)}-TrtrLG + ®[G] (16) 
which is stationary with respect to variations in G. Thus, 

8G' 

and the grand potential is constructed with <&[G] as the building block such that it 
is stationary for variations of the spinor G that satisfies the exact matrix Dyson 
equation, eq.(15). Here tr stands for the same notation as in Baym [14], namely, 

trAB= \ dl \d2A{l2)B[2l + ), (16a) 

o o 

/3 is the inverse temperature. Also, as in Baym, the choice of the branch of the 
logarithm is such that the variation of the first term in eq.(16) is taken to be of the 
form 8{tr\n(-G)} = -tr{(8G~ 1 )G] .This formulation of the many-body theory is 

known in the literature as the -derivable method and has recently been 
generalized to include the density functional formalism in [14a,b]. With this 
introduction of the one-particle spinor-Green function and the corresponding 
spinor- self-energy of the system under consideration, we are now in a position to 
generalize the Fukutome symmetry considerations of the HF solutions and 
subsequent classification of the broken symmetry solutions of a magnetic system. 

Following [13], the most general forms of the spinor-Green function and the 
spinor-self energy function with no spin-symmetry considerations are expressed in 
terms of the Pauli matrices: 

G(\2) = ±{g n (\2)+lg s (\2)}, (18) 
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E(12) = {<x b (12)+t-^(12)}, (19) 
We first note that the physical electron density and the physical vector spin density 
are respectively given by 

n(l) = (n(l)) = -*TrG(ll + ) = -ig n (ll + ). (20) 

and f(l) = (f(l)) = -/7>{TG(ll + )} = -i| s (ll + ). (21) 

The corresponding physical particle current vector and the physical spin current 
tensor respectively are given by 

/«=W=-i*-^)G(" , )} M . 

1 XT XT 

t =-^{( v >- v ') &(11,) U- 

Ad) = (a(d) = -^M^ - ) G " 1 r >U 

-iffv.-v.^di')}^,. 

Here we have expressed these quantities of interest in terms of the corresponding 
spin-scalar and spin-vector components of the full spinor-Green function, eq.(18). 
From eq.(9), we have 

d_ 

dt x 



i r , r r , 



2 mi 



>G ( (11') 



(24) 



= -\d{2)< 



-^{(f-/(l))G<(ll + )-G<(ll + )(f./(l))} = 
£ > (12)G < (21 + ) + G < (12)£ > (21 + ) 

> 

-Z < (12)G ) (21 + )-G ) (12)i: < (21 + ) ' 

From this, using eqs. (18, 19), we deduce the following general continuity 
equations relating the densities and currents given in eqs. (20, 21, 22, 23): 

a>(12)g<(21 + ) + ... 
+ ^>(12)-g(21 + )+... 
- t ,r r , f f^(12)^<(21 + )+. 



<§i) 



r r 



+ V 1 .y(l) = -Jrf(2> 



(25) 



+^>(12)xf<(21 + ) 



(26) 



The dots in the above denote other terms arising from the indicated manipulations 
required to express eq.(24) in terms of the particle- and spin-densities. Eq.(25) is 
the usual continuity equation for the particle density while eq.(26) is the 
corresponding one for the vector spin-density. The third term in the left hand side 
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is the torque term due to external field acting on the spin vector while the right 
hand side involves interaction contributions which are of two kinds. The first one 
is due to spin vector modified by the particle density while the second one is the 
cross product of the spin vector with the field due to interactions. The difference 
between the itinerant and the localized spin cases thus becomes clear. In the 
itinerant case, the divergence terms and the others on the right side contribute to 
give rise to spin wave dispersion while in the localized case, the divergence term is 
absent (no spin current) and only the cross product of spins contribute in the right 
hand side. We will return to this in the next section. 



(c) Fukutome's Classification of the Broken Symmetry Adapted Green 
functions and the Corresponding Self Energies 

Consider the Dyson equation in the absence of external fields obtained by dropping 
w in eq. (15). It transforms under any unitary transformation, U, to the form 



G _1 = G _1 - E' , where G _1 = U 



dt x 2m C ^ ^ 



(27) 



G' = UGU\ and E' = t/Et/ 1 ' = E[G']. 
The Hamitonian, H, given by eq.(l), leads to the self-energy as well as the first 
term in the right hand side of eq.(27). Now G _1 = Gq 1 because it is the 
noninteracting part; the self energy, on the other hand, reflects the contributions 
due to interactions, and is a therefore a functional of G, and has the form given in 
eq.(ll). Two cases arise as with the HF theory [1]: (a) G'=G, i.e., G commutes with 
U and thus is invariant under the transformation, U, and E' = E. In this case, U 
represents the self-consistent symmetry of the corresponding E -scheme (compare, 
HF scheme [1]); (b) G is not invariant under U but leads to G' and E' obeying the 
same form of the equation as G. Then U represents a broken symmetry. Fukutome 
[2] pointed out in the context of HF scheme that if Ul belongs to a subgroup of the 

group of transformations U, then G'HF and the corresponding self-consistent E HF , 
form the broken symmetry adapted solutions! From this it is clear that such a 
scheme holds for general self energy that appears in the equation for the Green 
function as in eq.(15), which we now explore in detail in this work. We will now 
present this generalized version of the Fukutome [2] analysis and describe these 
broken symmetry possibilities. 

We first consider only the group S of all spin rotations {j jJ }, I being 
the unit operator in the rotation group, and the group T of time reversal operation 
{l,T }which is of order two, I is the unit operator and 
T =iT£ , C = comp lex conjugation operator, under which the Hamiltonian (1) 
is invariant. As in [2], the inclusion of spatial symmetry has to be dealt with 
individually depending on the type of spatial symmetry one wants to consider. The 
general consideration involving only S and T is sufficient for the present and 
indeed, as will be discussed in Sec.IV, the results obtained here form the basis for 
the vector-spin-density-functional theory of magnetic systems, where the 
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interacting electron gas system provides the underlying functional for investigating 
the properties of many important inhomogeneous systems (see for example, [9, 
10]). Since S commutes with T, the group of all S and T is the direct product 
group, 3 = SxT = {l xljj JJ T }. There are eight different subgroups of 3 

which we will now enumerate. We note at once that IxJ , S, and T individually 
are three subgroups of G. The spin rotations A(e), about a fixed axis denoted by 
the unit vector e, is a subgroup of S: A(e)=p (e, 6)}, where p (e, 6)} is given 
by 

U («,0) = exp[-i0/2(f-e)] = cos|-/(f-e)sin|. (28) 

The fifth subgroup M( e' ) is of order two consisting of unit operator and the 
combined operation of T with a spin rotation through an angle K around an axis 
e' : M(e' )={l,TU (e' ,7i)}. There are two more subgroups that arise from the 
product of the elements A(e) and M(e'): 

A(e)xM(e')=p (e,0),U {e,6p {e\n)T } (29) 
This gives rise to two groups, one, when the unit vectors are orthogonal, 
e -e' = 0,and two, when e = ±e\ which is the group A(e)xT. Collecting them all 
together, we have the eight subgroups listed in Table I. 



TABLE I: THE EIGHT SUBGROUPS 

SxT S 
A(e)xT A(e)xM(e')(e-e' = 0) A(e) 

T M(e') Ixi 



In Table II, we have given the eight broken symmetry adopted forms of the Green 
function and the corresponding self-energy derived from the considerations given 
above. There is a corresponding set of self-energy functional that goes with these, 
which are also given in this Table. 
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TABLE II: BROKEN SYMMETRY ADAPTED GREEN FUNCTION 
SOLUTIONS OF THE DYSON EQUATION 



Invariance involving 
Time 

—> 

Invariance involving 
Spin 

i 


Group T of Time 
Reversal 


Group M(e') 
consisting of T and 71 

rotation about e' - 
axis. 


I 


Group S of all spin 
rotations 


(PARAMAGNETIC SYSTEM) 
G 1= l/2(g n ), with j = 0. 

Non-magnetic insulator. 


Charge current wave 
system 

G 2 =V2(g n ) 
with 7 ^0. 

Nonmagnetic metal or 
semiconductor 


Group A(e) of spin 

rotations about e- 
axis. (Axial) 


Axial spin current 
wave system 
G 3 = 1/2 

with 7=0 and 

z 3 = 

Insulating e — axis 
antiferromagnet 


Axial spin density 
wave system 
G 4 = 1/2 

with 7=0 and 
Is = 0. 

Insulating axial 
antiferromagnet. 

(e'.e = 0) 


Axial spin wave 

system 

G 5 = 1/2 

with 7^0 and 

z 5 = 

Itinerant electron axial 
antiferromagnet 


I 


General spin current 
wave system 

such that 
j = 0;s = 0,but 

fs*0- 

Insulating spin current 
density state 
(antiferro.) 


General spin density 
wave system 
G 7 = 1/2 

( gn+ (leyj + " 

Itinerant electron 
helical SDW state 
(e.g. Overhauser) 


General spin wave 
system (with no 
constraint on spin 
density vector) 

1 ^ IT 
G S = -(gn + *-gs)- 

Most general itinerant 
electron SDW state 
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In this Table, we have used the invariance of the continuity equations, eqs.(25, 26) 
under time-reversal operation which in turn lead to respective transformation 
properties of particle-, particle-current-, spin-, and spin-current- densities: 

n^n; j ->-/'; s -> -s; and j s -> j s . (30) 
This representation is equivalent to but more physical than the one given by 
Fukutome [2], who expresses his results in terms of the real and imaginary parts of 
the density matrix. It may be worth pointing out that several types of itinerant spin- 
density-wave states occurring in rare earth systems were discussed before in a 
phenomenological way [18]. 

III. LINEAR RESPONSE FUNCTIONS - SYMMETRY CONSIDERATIONS 



The well-known theory of linear response functions in itinerant magnetic 
systems may be expressed in terms of the functional derivatives of the inverse 
Green function, eq.(15), with respect to the external fields f a , of eqs. (6,7). From 
these one examines the possible collective excitations in the system, (see for 
example, [13,15]). Additionally these response functions allow us to test the 
stability of the state (for example, in [13, 15]) that is obtained as the solution of 
eq.(9 or equivalently 15). Alternately, the stability question may be reformulated as 
the criterion of minimum free energy, as was originally done in [14] and in a more 
physical way recently in [16] where a localized site-representation was used. Such 
a version is not applicable to itinerant electron system, where one has contribution 
from the mobile electrons as is evident from the work of [13,15]. One way to 
realize this is to note that in eq.(26) there is a contribution from the divergence of 
the spin-current which would be absent when there are no itinerant electrons. To 
bring out this feature, instead of repeating the work in [13, 15], we give here a brief 
outline of the linear response theory. The various types of response functions, % ap , 

can be expressed as appropriate variational derivatives, which in turn are expressed 
in terms of the corresponding appropriate variational derivatives of the self-energy: 

8s (I) 8g(u + ) 



8f{2) a 8f{2) 



-iTrtr\r a GA p G; A,(12;3) = = S(12)S(13)r p + J^g. 

Recalling the definitions given in eqs. (7, 20, 21), we have 16 response functions, 
involving the particle density and the three components of the spin vector density. 
In the localized electron scheme, the cross terms involving the density and spin 
vector will not appear. This is another important difference between the localized 
and itinerant electron systems. From such a linear response theory, one often 
deduces the low-energy collective excitations in the system, such as spin waves in 
ferro-, antiferro- and SDW- systems. Such a discussion may be found, for example, 
in [15]. 
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The objective in [16] was to derive expressions for the energy of the known 
localized effective spin models of magnetic systems. From this one can also 
deduce the low energy collective excitations in the system. In the itinerant electron 
system, it should be pointed out that there are important contributions not found in 
the localized spin systems, as is evident from [15], for example. We adopt the 
alternate formulation in this section to obtain a generalization of the "local force 
theorem" in [16] which is applicable to both the localized and the itinerant systems. 

Unlike in the general linear response theory described above, we now 
consider the effect of an infinitesimal rotation of the spin about a general direction 

denoted by a unit vector e, 86 = £80, obtained from eq.(28): 
Su = (l-iS0-z/2), on the free energy, eq.(16), in the same fashion as in [16] 
holding G fixed: 



8Q = -TrtrG(8Z); 

r r 
^ 80 , r _ r\ 80 ( r , r r x , r r x r\ 
8L = i — -(zZ-Zz) = i — -(z(z-g s )-(z-cj s )z) 



= - T. 1 80 x <j s ), and so 



(31) 



XT XT XT 

8Q = 8dtr(S s xg s ) = 8dV. 
Here we used the representations given in eqs.(18, 19) and some well-known 
identities to arrive at the expression for the torque, V, due to the rotation. This 
expression differs from the one given in [16] because of a factor of half in our 
definition of the spinor-Green function in contrast to theirs. As in [16], the Dyson 
equation gives us the sum rules on the components of the spinor Green function, 
after using eq.(15): 

GG" 1 = 1 = G^Gq 1 - L) = Ug n + I ^(Go 1 -o n -z- a s ); 



G~ l G = 1 = (Go 1 - L)G = (Go 1 - a n - I S s )^{g n + I g s ). 

From which we get the relations: 

( -\ \ r r 
[ G o -°n)gn- (J s-gs = 2 - 

/ \ XT XT XT XT 

(Go 1 - o n )8s - °s8n ~ i(J s x 8s = 0- 



(32) 



(33a,b) 



(34a,b) 



8n{G Q l -<y n )-g s -o s = 2. 

8s{ g q 1 -<* n )- g n o s - ig s xg s = 0. 

From these we derive the following sum rules corresponding to those given in [16] 
but in a coordinate representation valid for both itinerant and localized electron 
systems. Multiplying (33a) by g n on the left and (34a) on the right we obtain 
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2 §n =8n( G O l -< J n)§n-§n{S S -8s) 

r r ( 35 > 

= 8n( G 1 - a n)gn ~ {s S ' °s)s n - 

Multiplying (33a) by g s on the left and using (34b) we obtain 
2 8s = 8n°s8 n -8s{°s-8s) + i{8s x °s)s n 

= 8n°s8n ~ {§S ■ °s)Ss + ^ni^S X 8s)' 

The second expression is obtained by multiplying (34a) by g s on the right and 
using (33b). From eqs.(31, 36), we obtain the general expression for the torque 
vector 

V = rr(cr s x fo ) = --^{o- 5 x {-g n G s g n + (g s ■ g s -)g s - ig n (o s x g s ))\ (37) 

By integrating back the expression given by eq.(31), effectively the spin only part 
of the free energy may be rewritten as 
£l sp =trTr(GZ) = tr(g s G s ) 

(38) 

= -^ tr {{ss-^s){8s-^s)-{^s8n)i^s8n)-i{^s x 8n^s)-8s)\- 

Here eq.(36) was used in further simplification. The second approximation symbol 
is because we have dropped the spin-independent contribution, trg n a n , arising in 
the first expression. Note however, that from eq.(35), g n ,o n depend on the spin 
vector. Had we kept this contribution, we would have additional contributions due 
to particle density-density and particle density and spin density vector mentioned 
in the linear response theory at the beginning of this Section. For completeness, we 

give this here:£2„, = trg n o n = ^{g^Gq 1 - G n )g n - o n g n (o s -f^J.Here eq.(35) 

was used in further simplification. It should be noted that the notation tr used here 
is as defined before in eq.(16a). In the site-local representation used in [16], eq.(38) 
goes over to that given there. 

Using Table II, we have eight types of Green functions and their 
corresponding self-energies associated with the allowed broken symmetry 
solutions. Using these in the general expressions for the spin-spin interaction 
energies obtained above, one can deduce the structures of the corresponding spin- 
spin contribution to the free energy of the system. 

IV. IMPLICATION TO VECTOR - SPIN - DENSITY - FUNCTIONAL 
THEORY 



The development of the vector- spin-density-functional theory especially the 
local (LSD) approximation, has in recent years produced a much better theoretical 
understanding of itinerant magnetism. (See for example [5,9, 10].) In particular, in 
[15], a spinor-Green function version of the LSD theory was developed where the 
self-energy in eq.(9 or 15) is taken to be of the form, local in space and time, 
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Wl2)s{v JHr [n,j;r 1 ] + T.W w [ii,s;r 1 ]}5(12), (39) 
where V xc is the spin-scalar part of the self-energy, which is in general a local 
functional of particle density, n, and spin density vector, s, whereas W xc is its 
spin-vector counterpart, for describing the itinerant magnetic systems. The 
approximate forms for these functionals arise as the respective functional 
derivatives of the exchange-correlation energy, E xc [n,s ] = E xc [s a ], (using the 
notation of eq.(7)) of a spin-polarized homogeneous electron gas (see references in 
[5], for example). In fact, V xc = 8E xc /8n, W xc = 8E xc /8s, or more generically 

written as V x ^ = 8E xc /8s a , so that JL LSD = T a V x ^ The resulting equation is the 
LSD equation, 



2 



r 



r 



i^-+^-v c (i)- Wn (i)-i-f(i) 

dt x 2m 



G(12)- 

(40) 



- {V xc [nJ;l] + f • W xc [nJ;l]}G(l2) = 8(12). 

This is solved self-consistently by numerical methods. 

The continuity equations for the particle- and spin-density derived from this 
LSD eq.(40) now take the form 

r (41) 

c ^+v l js(i) = w xc xi(\). 

The second equation for the spin-density vector has a contribution from the torque 
due to the spin polarization of the spin system, besides having a contribution from 
the divergence of the moving spins producing a spin-current density. 
Proceeding as in eqs.(33, 34), we have the equations 

(Go 1 -V xc )g n -(w xc )-g s = 2. 



{Go 1 ~ V xc )g s - (W xc )g n - i(w xc )x l s = 0. 

g n (Go l -V xc )-g s -(w xc ) = 2. 

8s{Go l -V xc )- g n (w xc )- ig s x (W xc ) = 0. 
From these we obtain the expressions 

2 8n = 8n( G 1 - V xc )g n ~ g n (W xc ■ g s ) 
= Sn{ G Q 1 - V xc )g n ~ (g s ■ W xc )g n . 



(42a, b) 



(43a,b) 



(44a,b) 
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2 8s = 8 n (w xc )g n + i(g s x W xc )g n - g s (w xc ■ g s ), 

, r r r , , r r , r (45a,b) 

= 8n ( W xc )8n + I8n ( W xc X 8 S J " [8 S ' W xc J<?S • 

Finally the expression for the spin-only contribution to the free energy in the LSD 
case is obtained in a manner similar to the one given in Sec.III, eq.(38): 

a T ='\ tr {(ss ■ W xc ){8s -W xc )- {8 n W xc ) ■ [g n W xc ) - ig n {W xc xg s \ W xc }. 

(46) 

The corresponding particle density and spin density contribution to free energy is 
"Jf = ^« = 2 fr { V ^( G o 1 -^k B -^ B (^-| s )}. In eq.(46), as in 

eq.(38), the first term represents a multi-spin interaction energy involving four 
spins or more, the middle term is like the spin-spin interaction involving two spins 
or more, whereas the last term is a Dzialoshinskii-Moriya - type interaction, 
involving three or more spins. In [15], a spatially slowly varying approximation 
was considered in an approximate way and only the middle term resembling an 
effective Heisenberg spin-spin interaction energy was derived along with the cross 
terms containing particle density and spin density. The higher order spin terms and 
the Dzialoshinskii-Moriya terms did not appear in that derivation because the self- 
consistent relations as in eqs.(44, 45) were not invoked. The above remark arises 
from the observation that the vector part of the self-energy functional, W xc , is an 
odd functional of the spin vector beginning with a linear functional of the spin 
vector. 

In Table III given below, the eight types of Green functions and the 
corresponding self-energies are given associated with the various broken symmetry 
types for the case of LSD theory. From the above general expressions for the free 
energy, one may then deduce the structure spin dependent energies that follow for 
each of these cases. The phenomenological description given [18] is here 
supported from the considerations of group theory. 

Applying the linear response theory outlined in Sec.III to the LSD scheme, 

8Z SL Ss Y 8 2 E 8s y 

one obtains a tensor: = tr — -^ {chain rule) = T a tr &g ^ y ^ . as was 

shown in [15]. The broken symmetry considerations leading to Table III may be 
applied to this tensor to deduce the corresponding eight structures. As shown in 
[15] and more recently in [10], the use of homogeneous electron gas results in LSD 
requires a subtle and important modification in incorporating the vector nature of 
the spin density, when studying the spin wave properties of itinerant magnets. 
Another way of expressing this point is that the traditional electron gas theory 
leads to Ising-like treatment of the spins which is converted into a Heisenberg-like 
treatment (see for example, [9]) by a spin rotation. This does not lead to correct 
answers as was shown in [10] and this is due to the subtle nature of the treatment 
of the spin vector in the theory. In [10], a perturbation theory approach has been 
presented to include this feature, thus making a significant difference. 
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TABLE III: BROKEN SYMMETRY ADAPTED GREEN FUNCTION 
SOLUTIONS OF THE LSD EQUATION 



Invariance involving 
Time 

— > 

Invariance involving 
Spin 

i 


Group T of Time 
Reversal 


Group M(e') 
consisting of T and 71 
rotation about e' - 
axis. 


I 


Group S of all spin 
rotations 


(PARAMAGNETIC SYSTEM) 
G 1= l/2(g n ), with J = 0. 

£1 = S(12)V XC 
Non-magnetic insulator. 


Charge current wave 
system 

G 2 =}/2(g n ) 

with 7^0. 
Z 2 = 8(12)V XC 
Nonmagnetic metal or 
semiconductor 


Group A(e) of spin 

rotations about e- 
axis. (Axial) 


Axial spin current 
wave system 
G 3 = 1/2 

(g n +(T-e)gi f ) 
with 7=0 and 

Z 3 = 5(12) 

(V xc + (le)w'J c ) 

Insulating e —axis 
antiferromagnet 


Axial spin density 
wave system 
G 4 = 1/2 

(gn+(*-e)g'J) 
with 7=0 and 

Is = 0. 
£ 4 = 5(12) 

{V xc + (x-e)w'j c ) 

Insulating axial 
antiferromagnet. 

(e'.e = 0) 


Axial spin wave 

system 

G 5 = 1/2 

_1_ 

with 7^0 and 

L 5 = <5(12) 

(v xc + (le)W f J c ) 

Itinerant electron axial 
antiferromagnet 


I 


General spin current 
wave system 

1 IT 

G 6 =-{g n + T.g s ) 
such that 

j = 0;s = 0, but 

t 

Js * 0- 
Z 6 =<5(12) 

(V XC + T.W XC ) 

Insulating spin current 
density state 
(antiferro.) 


General spin density 
wave system 
G 7 = 1/2 

( gn+ (leyj + " 

Z 7 = (5(12) 

f V xc+ (le<)w'' c+ \ 

Itinerant electron 
helical SDW state 
(e.g. Overhauser) 


General spin wave 
system (with no 
constraint on spin 
density vector) 

1 ¥ IT 
G 8 = -{gn + T-gs)- 

22 8 = 5(12) 

(V XC + T.W XC ) 

Most general itinerant 
electron SDW state 
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V. SUMMARY AND CONCLUDING REMARKS 



In summary, the structures of the Green functions and their associated self 
energy arising from group theoretical considerations of the spin rotation and time- 
reversal invariance are given in Table II for general magnetic many-electron 
systems and in Table III, for the vector-spin-density functional formalism of the 
itinerant electron systems. This is expected to systematize the procedure in the 
analysis of magnetic structures that may appear in magnetic nanometric systems 
and in magnetic atomic clusters, just as the earlier similar work of Fukutome did 
systematize the Hartree-Fock solutions of magnetic states of molecular systems. 
We also consider consequences of this by setting up the linear response theory and 
an alternate version of it in the form of "effective spin Hamiltonian", to exhibit the 
differences between localized electron systems with those where the electrons are 
itinerant. We hope that the search for various types of magnetic structures in 
nanometric [8] and atomic cluster [11] systems, particularly when the clusters are 
deposited on semiconductor substrates [12], will be systematized by the procedure 
given here. In these systems, as seen from the work in [1 1], magnetic properties are 
correlated with the geometric structures of the clusters; the group theoretical 
analysis presented here is expected to make the search for this feature systematic. 
In this context, the work presented here incorporated into that in [17], may be 
expected to lead to efficient procedures for computation of magnetic properties in 
nanometric systems and in atomic clusters. 

We have here pointed out the significance of incorporating the vector nature 
of spin density in LSD theory, particularly in the modification needed in the 
traditional use of the homogeneous electron gas results. In this context, the recent 
work in [10] should be mentioned as an important step in a proper treatment of the 
vector spin. In [10], the transverse part of the vector-spin was incorporated in a 
perturbative way and was shown to lead to a better understanding of magnetism of 
Iron than previously [9]. The localized treatment of spin interactions deduced in 
[16] is here generalized to itinerant electron magnetic systems and the differences 
arising from this are spelled out. It may also be pointed out that the 
phenomenology of the various types of SDW structures given in [18] may be 
deduced from Table III. 

In this paper, we have not included the gauge group needed to incorporate 
superconducting phases nor have we included the lattice translation group. The 
addition of the translation group into the considerations given here brings in the 
irreducible representation characterized by the q-vector associated with the 
Brillouin zone and the little group of the q-vector. The gauge group is also useful 
particularly because of the newly discovered high Tc superconductors involving d- 
and s- wave pairing possess interesting vortex structures. Inclusion of these within 
the mean field (HF) theory has been reported in some special cases [2, 4]. The 
generalization of all these features in the Green function framework will be 
addressed in a future communication. 
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